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■ We calculate the first moment of the photon structure function, (x) 7 = 

Jq 1 dxF^ix, Q 2 ), on the quenched lattices with /3 = 6.0 using the formalism 
developed by the authors recently. In this exploratory study, we take into 
£> \ account only the connected contractions. The result is compared with the 
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Abstract 



experimental data as well as model predictions. 



The parton (quark and gluon) distributions in hadrons, which can be measured from 
lepton-hadron deep inelastic scattering and other hard processes, have played a crucial role 
in understanding high-energy scattering and the hadron structure. Recently it has become 
possible to calculate moments of these distributions from first principles using lattice QCD 

There are many other physical observables which are yet to be calculated from lattice 
QCD. A special class of of these involve either virtual or real photons. Until recently, 
the conditions under which these quantities can be calculated using lattice QCD remained 
unclear. The reason is that the photon is not an eigenstate of QCD. Rather, a "photon" state 



in nature is a superposition of the U(l) gauge boson and quark-gluon configurations which 
are suppressed by the electromagnetic coupling. The matrix elements of QCD operators in 
the photon states are time-dependent correlations which are defined in the Minkowski space, 
and the standard method used for calculating hadron matrix elements on the lattice is not 
directly applicable |j]-f§. 

In our earlier paper |4[], we showed that the matrix element of a quark-gluon operator 
between photon states can be evaluated using lattice QCD. The relevant expression is 

( 7 (pA / )|C(0)| 7 (pA)) = -e 2 J d 4 xrfV" (a;4 - ^/4) e- ^^? • (^; - ^;) (0|T i^ e*(A , ) • J(x)O(0)e(\) ■ J(y)\0) 

(1) 

where every quantity has been expressed in the Euclidean space: Te is the Euclidean time- 
ordering; the electromagnetic current J M (// = 1,2,3,4) is defined as Y^f e }'4 > }l^' l l ) f with 
Euclidean 7 matrices 7, = —i^ l M {i = 1, 2, 3) and 74 = 7 ; the Euclidean photon polarization 
vector is defined as = e l M and e 4 = ie°(e* 4 = ie*°); uo is the Euclidean-space photon energy. 
[In Ref. 0, we did not continue the electromagnetic current and the polarization vector to 
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Euclidean space and superficially the result there differs by an overall sign.] The above 
expression is the same as the naive analytic continuation of the matrix element from the 
Minkowski space. However, Eq. ([!]) is only valid when there exists an energy gap between 
the photon energy and the lowest hadronic state of the same quantum number. Otherwise, 
the integral becomes divergent. 

One of the quantities we can study using the above formalism is the photon structure 
functions, which can be measured from collisions between real or virtual photons with highly 
virtual ones, achievable in e + e _ — ► e + e _ +hadrons fHH. While a photon is normally con- 
sidered a structureless particle, it can fluctuate into a charged fermion-antifermion pair or 
more complicated hadronic states which can be revealed through interactions with a highly 
virtual photon. In fact, the photon structure functions can be defined from its hadronic 
tensor and the moments of the structure functions can be obtained from operator product 
expansion, much in the same way as for the hadron structure functions. There exists a sub- 
stantial amount of experimental data for unpolarized structure functions already ||; and 
future experiments from HERA and e + e~ are expected to measure the polarized structure 
functions as well |7]]. (For a recent review on theoretical and experimental progress on this 
subject, refer to Ref. ||.) 

Here we report on the first lattice study of the structure function F^. More specifically, 
we measure the first moment of which corresponds to the average fractional momentum 
carried by the quarks (weighed with the fractional charge squared for each flavor) in a real 
photon. The measurement is done on /3 = 6.0 lattices, corresponding to Q 2 ~ aT 2 ~ 
(2.4 GeV) 2 , in the quenched approximation for three degenerate flavors of quarks. Also, we 
choose to ignore the contributions from disconnected diagrams in which the quark lines do 
not connect any of the external vertices. 

We begin by briefly reviewing the definition of the photon structure functions and their 
relation to the cross sections. We closely follow |§ for the definition of F^x^Q 2 ). For the 
deep inelastic e + 7 process in a frame where the incoming electron and the photon are 
collinear, we write 



The relevant kinetic variables are Q 2 = —(k — k') 2 , u — p ■ q, q = k — k', x = Q 2 /2(p ■ q), 
y = q ■ p/k ■ p. The photon structural information is contained in the following tensor, 



e ± {k)+ 1 {p)^e ± {k')+X{p x ) , 
k» = (E,0,0,E) , 
k ,fl = (E',0,E'sm6,E'cos9) , 
^ = (£ 7 ,0,0,-£ 7 ) . 



(2) 
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p-q 



V 



Here J M (x) = 'Ylif e f' l l ) f( x )lii' l Pf(. x ) i s the electromagnetic current and e/ is the fractional 
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charge of the quark of flavor /. The covariant normalization has been used for the photon 
state, (7(pA)|7(p / A')) = (2n) 3 2uj^(p - p>)5 x , x , 

Since we are interested in the unpolarized photon properties here, we have averaged over 
the photon polarization A in the above equations. [The polarized photon structure functions 
are interesting by themselves and can be studied in a similar approach 0.] 

Using these variables, the differential cross section is given by 



da 



dE'd cos Q 4 E^ 
or, in terms of x and y and scaling functions 
da 



2E 2 W^(u, Q 2 ) cos 2 - + W?{y, Q 2 ) sin 2 - 



\Q-na 2 ra EE 1 



dxdy 



l-y)F2{x,Q?)+xy 2 Fl{x,Q 



(4) 



(5) 



where the dimensionless scaling functions are defined as F^x^Q 2 ) = vW^iv-iQ 2 ) and 
F?(x,Q 2 ) = W?(v,Q 2 ). 

In the parton model, the photon scaling function F^ix) is related to the quark distri- 
butions qj(x) in the photon — the probability of finding a quark with momentum fraction 

x, 

F^x)=x.^e}(q}(x)+q}(x)) . 
/ 

Taking into account QCD radiative corrections, the moments of the photon structure func- 
tions can be written in the following factorized form, 



f x n - 2 F2(x,Q 2 ) dx = J2c 2nJ (%,g( 
Jo j V ^ 



Oi) ) (x- 1 )} , (n = 2,4,...,) 



where Cmj are the perturbative coefficient functions and are the nonperturbative 

moments of the parton (quark and gluon) distributions (/ here refers to gluons as well). For 
quarks, the moments can be calculated as the matrix elements of local operators, 

\ ^2(liP^)\0 {fll ... tln }j\jipX)) = 2(x n - 1 y f [p fll ■ ■ - Pfin - traces] , 
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n-l 



(6) 



where Dfj,=D^ — is the covariant derivative and {■ ■ •} denotes the symmetrization of 
indices. 

In the following discussion, we ignore the radiative corrections to the coefficient functions 
and take C 2n j = e 2 for quarks and for gluons. We then define (x) 7 = F^{x, Q 2 )dx, 
which can be calculated as the matrix elements of the following operator; 



( r t) ( , 



f 



x) 



(7) 
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In lattice QCD, we need to define the corresponding lattice operator in Euclidean space: 
0<g?\x) = & f (x)r° u (x,y)My)+^ f (y)r° u (y,x)M' 



^{x 
f,y 

rJUz, V) = \ b„(U v (x)6 y>x+i> - Ut{y)5 v ^) + v)] , (8) 
where the lattice spacing(a) is set to 1 and the lattice covariant derivatives are defined as 

ip(x) = -(Ufa)if>(x + p)- Ufa - jj)il>(x - /*)) , 
^(ar) 5 M = ]M{x + p)Ufa) - ip(x - flUfa - p,)) . 
The corresponding Euclidean matrix element to Eq. is 

\ ^(7(pX)\0^\ 7 (pX)) = 2i( X y P , EPvE , (9) 

A 

where the photon Euclidean four- momentum is defined as p^E = (jP,iw)- 
The photon matrix element in lattice QCD is 

(l(p\')\O, u (0)\l(pV) ilaU) = -e 2 4(A'MA) 

x J^e^-^e-^-y^OlTEJ^O^^J^lO) . (10) 



The Euclidean Green's function can be calculated as follows, 

*././'./" 



\T E .Ux)o^\o)My)\o) = ~ E w e ? 



-E4(2 



K) 3 



(Tr[ la S(x, 0)1^(0, z)S(z, y)^S{y, x))) g (11) 



[Tr[ lp S{y, 0)1^(0, z)S(z, x) la S(x, y))) + (Tr{ la S(x, z)T° u (z, 0)5(0, y) lp S(y, x)}) 



+(Tr[ypS{y, z)T°(z, 0)5(0, x) la S(x, y)]) t 



+ ... 



where we have shown explicitly the connected contraction with one trace and the ellipsis 
represents disconnected diagrams with more than one trace, and (■ ■ - ) g denotes averaging 
over gauge configurations. The propagator S(x,y) is solved from '^2 x D(y,x)S(x, z) = 8 y>z 
where 

D(x, y) = 5 x , y Y)U,(x)5 y , x+fl + (1 + 7 ^(y)<W , (12) 

is the scaled Wilson-Dirac operator. The Dirac and color indices are suppressed for read- 
ability and the traces are over those indices. 



In this exploratory study, we will ignore the disconnected contribution. For the connected 
contraction, we can show, using the charge conjugation property of the Wilson-Dirac oper- 
ator, that the average traces for the left and the right derivatives are the same and, hence, 
we only need to evaluate one of them. The Wilson-Dirac operator D(x, y) is invariant under 
the charge conjugation: 

V>(x) -> C^(x) T , $(x) -> -*l){x) T C-\ U^x) -> (U^x))* , (13) 

where the matrix C satisfies C^j^C = — and T denotes transpose operation. Using 
Eq. (|13"D, we can show 

S(x,y\U) T = C- 1 S(y,x\U*)C , (14) 

where we explicitly denote the dependence of the propagator on the gauge field. 

Now we can relate the traces which arise from the left derivative of T® u with those from 
the right derivative. For example, 

T±[7 a S(x,0\U)jMr{0)S(hy\U)>yf,S{y,x\U)] 

= Tr [Siy^Wf^SifiMUfUli^SixMUfll] (15) 
= -Tr [S(x,y\U*)^S{y,fi\U*)(U%(0) lv S{0,x\U*) la ] . 

Using the previous equation and averaging over gauge configurations, we get 

(Tr[ la S(x,0)T^{0,z)S{z,y)^S{y,x)}) g = (Tr[^S(y,z)T^(z,0)S(0,x) la S{x,y)}) g , 

(16) 

and similarly for x <-> y and Finally, we arrive at the following expression used for 

actual calculation: 

(j(pX)\O„ v (0Mp\))^ = - e 2 e;(A')e,(A)^e^-^e-^^^ e 4 / (2 K ) 3 x 

x ,y zj 
2 • (Tr[ la S(x, 0)1^(0, z)S(z, yhpS(y, x)] + Tr[ 7p S(y, 0)^(0, z)S(z, x) la S(x, y)]) g . 

The above matrix element is evaluated by summing over the propagators starting from 
point sources placed at the position of O instead of one of the electromagnetic vertices. The 
reason for this strategy, different from hadronic structure functions [0-0], is that, as evident 
from Eq. (|ll|), the electromagnetic currents at both source and sink should be summed over 
all the time slices with proper momentum dependence, in contrast to the hadron structure 
functions where the time coordinates of the hadron operators are fixed. [Because of the 
translational invariance in the time direction, one can still fix the time coordinate of one 
of the electromagnetic vertices instead of the operator observable. The answer shall not be 
much different if the time direction is sufficiently large.] 

In practical simulations, two sequential propagators are generated from each of the 12 = 
(4 x 3) point sources 

M a (p,u,y) = J2 x S(y,x)7 a S(x,0)e-& s+ ^ , 

M' (p, co, x) = V S(x, y)j p S(y, 0)e+^™ . (17) 
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The complementary propagator generated from the operator O is 

NZ,(x) = 52rZ,(0,z)S(z,x) , 

z 

which can be calculated from the point sources as N®J(x) = J2 Z 75^( x ' z )l5^u (, z > 0)- I n 
terms of the above propagators, 

(^)conn./"em = ~ — ~ e}e*e„(2/c) 3 



x 2 • 



(18) 



where a and (3 are fixed indices (without summation) and the polarization vector can be 
chosen as a linear one in the x direction, e p = (1,0,0,0). For 3 degenerate quarks, 5^/ e / 
= (— 1/3) 4 + (2/3) 4 + (— 1/3) 4 = 2/9. [Alternatively, one can generate doubly sequential 
propagators: 

M {2)aP {z) = S(z, y)lpS(y, x) la S(x, Q yv-y-^ e -^ + ^ ; (19) 

x,y 

M{ 2)aP {z) = Y,S(z,x) la S(x,y) lp S(y,0)e^'^e-^ + ^ , 



which would eliminate the need for separate inversions for different operators. However, the 
successive inversions are then not desirable numerically when high precision is needed.] 

We use (3 = 6.0, 16 3 x 32 quenched configurations generated by Kilcup et al. ||, available 
from NERSC lattice archive. The lattice spacing is determined by taking the p meson mass 



at the chiral limit: a ~ 2.4 GeV [|K|. To decide on the boundary condition for the Wilson 
fermion in time direction, we made a test run with p = 5.7, 8 2 x 16 x 32 lattices. For the 
antiperiodic boundary condition, the finite size effect gives a noisy (x) 7 at a m q ~ 200 MeV 
(k = 0.16). On the other hand, the effect is under much better control for a m q ~ 100 MeV 
(ft = 0.165) for the Dirichlet boundary condition. For heavier quark masses, both boundary 
conditions give consistent results. Therefore, we use a Dirichlet boundary condition to obtain 
the main result of the paper. 

The photon momentum fraction (x) 7 is measured by evaluating ( r y(p\)\O 3 4(0)\'y(p\)) 
for the momentum k = (0, 0, 27r/16), to = 27r/16. We use hopping parameters k = 153, 
0.154, 0.155, which corresponds to pion masses ~ 0.424, 0.364, 0.301 in lattice units, 
respectively |TD[. We also evaluate (x) 7 for k = (0,0,7r/16), uj = 7r/16 (k = 154, 0.155), on 



the same Monte Carlo lattices, which is possible by employing the antiperiodic boundary 
condition for the pseudofermion fields in the z direction. 

To be close to the real world, the moment is calculated for three flavors of quarks with the 
degenerate masses. In the quenched approximation, the "p meson" is the lowest hadronic 
state with photon quantum number. With these parameters, m p = 0.426 (in lattice units) 



for the lightest quark mass (k = 0.155) [[H]] and \lm 2 p + k 2 = 0.577. This gives L t x 



m 2 p + k 2 — uj 
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To test if this is large enough, we measured (x) 7 with 3 different time positions of the 
operator O. Figure 1 shows the value of the first moment of F^ for k = 0.153, 0.154 with Tq 
= 16,17,18 (time coordinate runs from to 31). The numerical values of (x) 7 for different 
To's are consistent with each other (Fig. |IJ), which shows the temporal size of the box 
is large enough for the mass range. Note that the lattice result has been converted to a 
continuum renormalization scheme (MS in this case) by 

0™(Q 2 ) = Zo ■ ^V) (20) 
Zo = 1 + 9 \ N \ T l (V^MQV) - (£ LATT - B m )) (21) 



16vr 2 2N c 

Here we use the renormalization constant calculated perturbatively for Q = a" 1 ||11|| , Zq = 
0.9892. (The mixing between the operator 0^ u (Eq. (^)) and gluon operator Ti[F^ p F piJ } is 
absent in the quenched approximation.) 

Figure ^ shows the first moment of K7 for 3 degenerate quark flavors. The values of 
(x) 7 /a em for both momenta (uo = 27r/16, 7r/16) are well within the error bar. A linear fit 
to the chiral limit (m^ = 0) gives ~ 0.72(8), which is significantly larger than existing 
experimental results 13. Also, existing theoretical models such as the GRV model [13|] and 



the Quark Parton Model (QPM) |TJ predicts (x) 7 /a em ~ (0.3-0.4) for the value of Q 2 
investigated here. 

We notice that the previous lattice studies of hadron structure functions [l|-|3| showed a 
consistent overestimation of the first moment of the structure functions {(x)). It has been 
suspected that the absence of sea quarks enhances the valence quark contribution, resulting 
in larger values of (x). However, recent results from unquenched structure function studies || 
show no significant deviation from quenched results, which suggests that much smaller quark 
masses and/or a larger lattice size, quenched or unquenched, may be necessary to approach 
the continuum value. So some amount of overestimation for (x) 7 is not unexpected. Still, the 
difference is more significant than those for hadron structure functions. This may suggest 
that F2 is actually larger at the extremes of x, where the experimental data is not available. 

Another source of discrepancy is that the contribution from disconnected diagrams may 
be large. We are currently investigating this. [It should be noted that for 3 degenerate 
(up, down and strange) quarks, the trace of electromagnetic operator, Tr[Jg(x)], is zero 
since Y2f e f = 0; we would only need to evaluate the diagrams similar to disconnected 
insertion (DI) diagrams studied for quantities, such as the strangeness magnetic moment of 
the nucleon JH|.] 

In summary, the formalism necessary to compute the moments of the photon structure 
function is presented. We then compute the first moment of the unpolarized photon structure 
function F^ on quenched (3 = 6.0 (Q ~ a -1 ~ 2.4 GeV) lattice configurations. The result 
is somewhat larger than the existing theoretical models and experimental results. The 
discrepancy is likely due to the quenched approximation and the disconnected contribution. 
If the result persists, it may indicate F^x) is significantly higher at near x = or x = 1 than 
the theoretical models show. However, the systematic errors present in the measurement 
need to be understood to make more meaningful comparison. 

This work is supported by funds provided by the U.S. Department of Energy (DOE) 
under grant No. DOE-FG02-94ER-40762. The numerical calculation reported here was 
performed on the Calico Alpha Linux Cluster at the Jefferson Laboratory, Virginia. 
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FIGURES 
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FIG. 1. The first moment of for /3=6.0 quenched configurations for the different time posi- 
tion of the operator (To). 
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m vs. <x> la 
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FIG. 2. The first moment of for [3=6.0 quenched configurations (a -1 ~ 2.4 GeV). 
to = tt/8 (vr/16) correspond to the periodic (antiperiodic) boundary condition for the pseudofermion 
fields in the z direction. m n is in lattice units (a = 1). 
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